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ABSTRACT 

We present KPNO 4-m and LBT/MODS spectroscopic observations of an H II region in the nearby 
dwarf irregular galaxy Leo P discovered recently in the Arecibo ALFALFA survey. In both observa- 
tions, we are able to accurately measure the temperature sensitive [O III] A4363 line and determine a 
"direct" oxygen abundance of 12 + log(0/H) = 7.17 ± 0.04. Thus, Leo P is an extremely metal defi- 
cient (XMD) galaxy, and, indeed, one of the most metal deficient star-forming galaxies ever observed. 
For its estimated luminosity, Leo P is consistent with the relationship between luminosity and oxygen 
abundance seen in nearby dwarf galaxies. Leo P shows normal a element abundance ratios (Ne/O, 
S/O, and Ar/O) when compared to other XMD galaxies, but elevated N/O, consistent with the "de- 
layed release" hypothesis for N/O abundances. We derive a helium mass fraction of 0.2509^o'oi23 
which compares well with the WMAP + BBN prediction of 0.2483 ± 0.0002 for the primordial helium 
abundance. We suggest that surveys of very low mass galaxies compete well with emission line galaxy 
surveys for finding XMD galaxies. It is possible that XMD galaxies may be divided into two classes: 
the relatively rare XMD emission line galaxies which are associated with starbursts triggered by infall 
of low-metallicity gas and the more common, relatively quiescent XMD galaxies like Leo P, with very 
low chemical abundances due to their intrinsically small masses. 
Subject headings: galaxies: abundances - galaxies: dwarf - galaxies: evolution 
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Metallicity determinations of the interstellar media in 
galaxies allow us to assess the chemical evolutionary sta- 
tus of a galaxy. In this regard, studies of the lowest 
metallicity, least chemically evolved galaxies are of spe- 
cial interest (Mateo 1998; Kunth & Ostlin 2000; Tolstoy 
et al. 2009; McConnachie 2012). Specifically, they are 
important for studies of: (1) star formation at low metal- 
licity, which is critical to understanding star formation in 
the early universe; (2) stellar properties of low metallic- 
ity massive and intermediate mass stars, which provide 
critical constraints on theories of stellar evolution; (3) 
constraints on the early enrichment of the pre-galactic 
medium; and (4) studies of the primordial elemental 
abundances. 

Giovanelli et al. (2013) have reported the discovery of 
a nearby dwarf irregular galaxy (AGC 208583 = Leo P) 
as part of the Arecibo Legacy Fast ALFA Survey (AL- 
FALFA, Giovanelli et al. 2005; Haynes et al. 2011). The 
ALFALFA survey is a blind survey in the HI 21cm line 
covering 7000 square degrees of high Galactic latitude 
sky. Leo P was discovered as part of a program to identify 
mini-halo candidates in the ALFALFA survey (Giovanelli 
et al. 2010; Adams et al. 2013). Follow-up broadband op- 
tical imaging by Rhode et al. (2013) clearly demonstrated 
the presence of a resolved stellar population at the loca- 
tion of the ALFALFA HI detection, and Ha imaging in 
the same study revealed a strong Ha source associated 
with a bright central star in Leo P. 

Here we report on optical spectroscopy of the emis- 
sion line source in Leo P. We present results of optical 
spectroscopy from the Kitt Peak National Observatory 
(KPNO) 4-m and from the Large Binocular Telescope 
(LBT) in Section 2, and derive H II region physical con- 
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ditions from the spectra in Section 3. We derive chemical 
abundances from each of the two spectra in Section 4. 
We then compare the results of these abundance anal- 
yses to similar analyses of low metallicity H II regions 
in Section 5, and discuss the implications for our under- 
standing of metal poor galaxies in Section 6. Finally, we 
summarize our conclusions in Section 7. 
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Fig. 1. — The spectrum of the H II region in Leo P taken with the 
KPNO 4-m. Note the clear detection of the temperature sensitive 
[O III] A4363 emission line and the very weak strength of the [N II] 
AA6548,6584 lines, indicative of a very low abundance. Ha, H/3, 
and O III A5007 are all off scale to allow better visibility of the 
weaker emission lines. 



2. OPTICAL SPECTROSCOPY OF LEO P 
2.1. The KPNO 4-m Spectrum 

The H II region discovered by Rhode et al. (2013) has 
coordinates of R.A. = 10:21:45.1 and dec. = +18:05:17.2 
(J2000). A spectrum of the Leo P H II region was ob- 
tained with the KPNO Mayall 4-m on 2012 April 20 un- 
der relatively clear skies and approximately arcsecond 
seeing. The Ritchey- Chretien Focus Spectrograph was 
used with a Tektronix 2048 2 pixel detector (T2KA). A 
316 line mm" 1 grating (KPC-10A), 1.5" slit, and WG345 
blocking filter were used. This spectroscopic configura- 
tion yielded a dispersion of 2.78 A per pixel, a full width 
at half maximum resolution of < 6 A, and a wavelength 
coverage of 3500-8300 A. Bias frames, flat-field lamp im- 
ages, and sky flats were obtained. Combined helium, 
argon, and neon arc lamps were acquired for accurate 
wavelength calibration. The spectrophotometric stan- 
dard star Feige 34 (Massey et al. 1988) was observed be- 
fore and after the observations of Leo P to yield a secure 
flux calibration. 

Four 1200 second exposures of Leo P were obtained 
with the slit at a fixed position angle which approximated 
the parallactic angle at the midpoint of the observations. 
This, in addition to observing at airmasses of less than 
~1.2, served to minimize the wavelength-dependent light 
loss due to differential refraction (Filippenko 1982). 

Standard procedures within the IRAF 10 ccdred and 

10 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



specred packages were used to bias-subtract, flat-field, 
and illumination-correct the raw data frames. The 
lacos_spec routine (van Dokkum 2001) was used to 
remove cosmic rays from the images; this process was 
checked carefully to ensure that sharp emission lines in 
our nebular spectra were not accidentally misidentified 
as cosmic rays by the software. The spectra were then 
extracted from the 2D images using the IRAF apextract 
package. The individual ID spectra were wavelength cal- 
ibrated by applying the solution from the HeArNe com- 
parison lamps, and flux-calibrated using the sensitivity 
curve derived from the standard star observations. The 
latter step also involved making airmass-dependent at- 
mospheric extinction and reddening corrections using the 
standard Kitt Peak extinction curve (Crawford & Barnes 
1970). The multiple fully-processed spectra were then 
combined into a single composite spectrum. 

Figure 1 shows the final calibrated 4-m spectrum. Note 
the clear detection of the temperature sensitive [O III] 
A4363 emission line which allows for a "direct" abun- 
dance determination and the very weak strength of the 
[N II] AA6548,6584 lines, indicative of a very low metal- 
licity (Denicolo et al. 2002). 
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Fig. 2.— The blue portion of the LBT/MODS spectrum of the 
H II region in Leo P. Note the clear detection of the temperature 
sensitive [O III] A4363 emission line and the detection of the He I 
A5015 line. The higher numbered Balmer series lines are labeled 
H16 - H7. Note that underlying absorption is not detected for the 
Balmer lines. O III A5007 is off scale to allow better visibility of 
the weaker emission lines. 



2.2. The LBT/MODS Spectrum 

After the KPNO 4-m spectrum revealed the very low 
oxygen abundance of Leo P, we re-observed the H II re- 
gion in Leo P in order to obtain a higher signal-to-noise 
spectrum and to detect weak emission lines over a larger 
wavelength range. A spectrum of the H II region in Leo P 
was obtained with the Multi-Object Double Spectro- 
graph on the Large Binocular Telescope (LBT/MODS, 
Pogge et al. 2010) on 2012 April 29 under relatively clear 
skies and approximately 0.6" seeing. The LBT/MODS is 
a double spectrograph with a dichroic that splits the light 
at ^5650 A and sends it to two separate spectrographs. 
LBT/MODS1 was used in longslit mode with a 1.0" slit 
imaged onto two 3072 x 8192 format e2v CCDs with 
15^m (0.12") pixels. The blue side of MODS1 covers a 



3 



wavelength range of 3200 to 5650 A with a 400 1 mm 
grating providing a resolution of 2.4 A. The red side of 
MODS1 covers a wavelength range of 5650 to 10000 A 
with a 250 1 mm -1 grating providing a resolution of 3.4 
A. Bias frames, flat-field lamp images, and sky flats were 
obtained. We observed BD+33°2642 and Feige 34 as the 
standard stars with a 5x60" spectrophotomctric slit mask 
near the parallactic angle. These are from Oke (1990), 
using the HST CALSPEC flux tables (e.g., Bohlin 2010) 
which extend further into the UV and Near-IR than the 
original Oke fluxes. Feige 34 was observed at an airmass 
that more closely matched the airmass of the target ob- 
servations, so it was used as a single calibration source 
instead of being used in combination with BD+33°2642. 
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Fig. 3.— The red portion of the LBT/MODS spectrum of the 
H II region in Leo P. Note the very weak strength of the [N II] 
AA6548,6584 lines, indicative of a very low abundance. Also note 
the detection of the [S III] lines at AA9069, 9532 and the emission 
lines from the Paschen sequence which are labeled P14 - P8. The 
noise spikes at wavelengths above A8500 are due to residuals from 
subtraction of the atmospheric emission lines. Ho is off scale to 
allow better visibility of the weaker emission lines. 

Three 900 second exposures were taken at an airmass 
values less than <~1.1, with the slit at a fixed position an- 
gle which approximated the parallactic angle at the mid- 
point of the observation. Wavelength calibrations were 
derived from Hg(Ar), Ne, Xe, and Kr lamp spectra taken 
with the telescope at zenith. MODS is equipped with a 
closed-loop active flexure compensation system, so only a 
small (< lA) residual flexure correction, calculated from 
night-sky lines in the spectra, was applied to the final 
wavelength calibration. 

Standard procedures within the IRAF 2-D package 
(twodspec) were used to bias-subtract, flat-field, and 
illumination-correct the raw data frames. Airmass de- 
pendent atmospheric extinction was corrected using the 
Kitt Peak (Crawford & Barnes 1970) and Cerro Paranal 
(Patat et al. 2011) extinction curves corrected approxi- 
mately to the 3220-meter elevation of Mt. Graham adopt- 
ing an atmospheric scale height of 7 km. 

The multiple sub-exposures were combined, eliminat- 
ing cosmic rays in the process. The resulting images were 
then flux-calibrated using the sensitivity curve derived 
from the standard star observation. Finally, a trace fit 
to the continuum source in the slit was used to extract 
the galaxy emission within an aperture that encompassed 



> 99% of the light. 

In Figure 2 we present the calibrated blue spectrum 
and in Figure 3 the calibrated red spectrum. Figure 2 
shows the clear detection of the Balmcr emission line se- 
quence down to H12 at A3750 due to the relatively high 
sensitivity of LBT/MODS at wavelengths below A4000 
(and continuing down to the atmospheric cut-off) and 
the relatively high spectral resolution. Note also the high 
significance of the temperature sensitive [O III] A4363 
emission line, the detection of the He I A5015 line, and 
the lack of obvious underlying absorption for the B aimer 
lines. Figure 3 shows the very weak strength of the [N II] 
AA6548,6584 lines, indicative of a very low abundance, 
detection of He I A6678 at high significance, and the de- 
tection of the [S III] lines at AA9060, 9532 with emission 
lines from the Paschen sequence. 

3. EMISSION FLUXES AND PHYSICAL CONDITIONS 
3.1. Emission Line Measurements 

Emission line strengths were measured from the two 
spectra using standard methods available within IRAF. 
In particular, the splot routine was used to analyze the 
extracted one-dimensional spectra and to measure emis- 
sion line fluxes by integrating under the line. In cases 
where emission lines were blended, multiple Gaussian 
profiles were fitted to derive the integrated line fluxes. 
Special attention was paid to the Balmer lines, which, in 
some extragalactic H II regions, are located in troughs of 
significant underlying stellar absorption (e.g., Berg et al. 
2011). For Leo P there is very little evidence of under- 
lying absorption. For the bluest Balmer lines, multiple 
component fits were attempted in which the absorption 
was fit by a broad, negative Lorentzian profile and the 
emission was fit by a narrow, positive Gaussian profile. 
However, even in these cases, the correction for underly- 
ing absorption was small, i.e., comparable to the uncer- 
tainty in the flux in the line, so simple integration under 
the line was used. 

The errors of the flux measurements were approxi- 
mated using 

<Jx « \j(2 x x 2 + (0.02 x -F\) 2 , (1) 

where N is the number of pixels spanning the Gaussian 
profile fit to the narrow emission lines and £ is the rms 
noise in the continuum determined as the average of the 
rms on each side of an emission line. For weak lines, 
the uncertainty is dominated by error from the contin- 
uum subtraction, so the rms noise term determines the 
uncertainty. For the lines with flux measurements much 
stronger than the rms noise of the continuum (usually the 
Ha lines and often the [O III] AA4959,5007 doublet) the 
error is dominated by flux calibration and de-reddening 
uncertainties. In this case, a minimum uncertainty of 2% 
was assumed (based on the uncertainties in the standard 
star measurements, Oke 1990), and the right hand term 
above dominates the uncertainty estimate. 

3.2. Reddening Corrections 

The relative intensities of the Balmer lines are used 
to solve for the reddening using the reddening law of 
Cardelli et al. (1989), assuming A v = 3.1 E(B - V) 
and the theoretical case B values from Hummer & Storey 
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(1987) interpolated to the temperatures derived from the 
[O III] lines. We used a minimized \ 2 approach to solve 
simultaneously for the reddening and underlying absorp- 
tion based on the Ha/H/3, Hy/H/3, and HJ/H/3 ratios 
(cf. Olive & Skillman 2001). The reddenings derived 
for the two spectra are both very low with C(H/3) < 
0.1. The underlying Balmer absorption is consistent with 
zero for both spectra. The reddening correction for the 
LBT/MODS spectrum can be tested by comparing the 
corrected higher numbered Balmer lines to their theoreti- 
cal values. The H9 A3935 and H10 A3798 corrected fluxes 
are consistent with their theoretical ratios to H/3 of 0.074 
and 0.054, respectively. The derived values are listed in 
Table 1. Note that in Section 4.4 we will re-derive the 
reddening (and other parameters) from the combination 
of H and He lines. 

The foreground reddening in this direction is estimated 
to be Ay = 0.07 magnitudes (Schlafly & Finkbciner 
2011), which is very weak, corresponding to the high 
Galactic latitude of Leo P (+54°). The total reddening in 
the LBT/MODS spectrum of C(H/3) = 0.09 corresponds 
to an Ay w 0.2. This is slightly larger than the expected 
foreground reddening, indicating only a small amount of 
reddening internal to the galaxy, in accordance with its 
very low metallicity (e.g., Cannon et al. 2002, found a 
value of Ay = 0.13 magnitudes of internal extinction in 
I Zw 18). The relative flux estimates and corresponding 
errors for the two spectra, corrected for reddening, are 
listed in Table 1. 

Overall, there is good agreement between the rela- 
tive fluxes in the emission lines. A priori, some small 
differences are expected due to the difference in slit 
widths and/or a slight difference in slit placement be- 
tween the two observations. One small difference occurs 
for the fainter Balmer lines Hll (A3771), H10 (A3798), 
H9 (A3825), H8+Hel (A3889), and H7+[Ne III]. For these 
five emission lines, the KPNO 4-m fluxes are systemati- 
cally less (though not significantly) than those from the 
LBT/MODS spectrum. We believe this to be an effect of 
the very weak underlying H absorption (since observed H 
absorption is nominally constant in terms of equivalent 
width, the weakest, i.e., bluest, lines are the most sensi- 
tive to the presence of underlying absorption). Because 
the resolution of the KPNO 4-m spectrum is roughly 
half that of the LBT/MODS spectrum, the broader un- 
derlying absorption will have a greater impact on the H 
emission line fluxes. The overall consistency between the 
abundances derived from the two spectra in Section 4 is 
a direct result of the good agreement between the two 
spectra. 

3.3. Electron Temperature and Density Determinations 

For the purpose of deriving nebular abundances, we 
adopt a simple two zone approximation, where t 2 and ts 
are the electron temperatures (in units of 10 4 K) in the 
low and high ionization zones respectively. For the high 
ionization zone, the [O III] I(AA4959,5007)/I(A4363) ra- 
tio was used to derive a temperature. We use the 5-level 
atom calculations available with the IRAF task temden. 
The derived temperatures are given in Table 2. We de- 
rive relatively high temperatures of 17,150ti3go K from 
the KPNO 4m spectrum and 17,350 K from the 

LBT/MODS spectrum, as expected for a low metallicity 



H II region. The agreement between the two independent 
measurements indicates that the H II region in Leo P is 
clearly very low metallicity. 

Because neither the [O'll] AA7320,7320 lines nor the 
[N II] A5755 line were detected at high significance, we 
cannot derive a temperature for the low ionization zone 
directly and, therefore, need to assume a temperature. 
We used the relation between ti and £3 proposed by Pagcl 
et al. (1992), based on the photoionization modeling of 
Stasihska (1990) to determine the low ionization zone 
temperature: 

t 2 ~ 1 = 0.5(i 3 _1 +0.8). (2) 

The [S III] A6312 emission line was only detected at the 
2 a level, so for the temperature in the [S III] zone we 
assumed a temperature based on the relationship derived 
by Garnett (1992): 

t(S +2 ) = 0.83(i 3 ) + 0.17. (3) 

The low and high ionization region temperatures are tab- 
ulated in Table 2. 

[S II] AA6717,6731 and [O II] AA3726,3729 were used 
to determine electron densities. These densities are both 
consistent with the low density limit. The [O II] mea- 
surement, which is allowed by the higher resolution of 
the LBT/MODS spectrograph, provides a stronger con- 
straint on the density due to the lower critical densities 
of the [O II] emission lines (Osterbrock & Ferland 2006). 
For all abundance calculations we assume n c = 10 2 cm~ 3 
(which is consistent with the 1 a upper bounds and pro- 
duces identical results for all lower values of n c ). 

4. NEBULAR ABUNDANCE ANALYSIS 

We analyze both spectra separately, compare the de- 
rived abundances, and provide adopted abundances and 
uncertainties in Table 2. Calculated errors in the abun- 
dances provide a statistical estimate only. Additional er- 
rors may be important, such as systematic errors due to 
temperature fluctuations (cf., Pena-Guerrero et al. 2012). 
However, adopting this approach means that the derived 
abundances and uncertainties will be directly compara- 
ble to most of those reported in the literature. 

4.1. Oxygen Abundance Determination 

We determine oxygen abundances based on our esti- 
mated two zone electron temperatures. Ionic abundances 
were calculated with: 

Mg) = lm m. (a) 

N(H+) I Hfj Jx{€) • 1 ; 

The emissivity coefficients, which are functions of both 
temperature and density, were determined using the 
IONIC routine in IRAF with atomic data updated as re- 
ported in Bresolin et al. (2009). This routine applies the 
5-level atom approximation, assuming the appropriate 
ionization zone electron temperature. The total oxygen 
abundance, O/H, is the sum of O+/H+ and 0++/H+. 

The oxygen abundance determinations for Leo P are 
given in Table 2. We derive an oxygen abundance of 12 + 
log(0/H) = 7.17 ± 0.07 from the KPNO 4-m spectrum 
and 7.17 ± 0.05 from the LBT/MODS spectrum. The 
results from the two spectra are in good agreement, and 
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we adopt the error weighted average of 12 + log(0/H) 
= 7.17 ± 0.04 as the measurement of the oxygen abun- 
dance. 

The oxygen abundance in Leo P is one of the low- 
est oxygen abundances ever derived for an H II region. 
Formally, the oxygen abundance in Leo P is equivalent 
to that of the famous low metallicity galaxy I Zw 18 
(with 12 + log(0/H) = 7.17 ± 0.04; Skillman & Kenni- 
cutt 1993; Izotov & Thuan 1999) and lower than that of 
SBS 0335-052E (log(0/H) = 7.33 ± 0.01; Izotov et al. 
1997). It is slightly higher than observed in SBS 0335- 
052W (log(0/H) = 7.12 ± 0.03; Izotov et al. 2005) and 
DDO 68 (log(0/H) = 7.14 ± 0.03; Pustilnik et al. 2005; 
Izotov & Thuan 2007) , which are the current record hold- 
ers. 

Note that there are some reports of even lower oxygen 
abundances in the literature, but these are typically from 
spectra with lower signal/noise ratios. A small overesti- 
mate of the [O III] A4363 line can produce artificially 
high values of the electron temperature, which translate 
into lower oxygen abundances. These objects are usually 
associated with derived electron temperatures well in ex- 
cess of 20,000 K, which is difficult to achieve in real H II 
regions, even with only trace amounts of oxygen present. 

The significance of the low oxygen abundance in Leo P 
will be discussed further in Section 6. 

4.2. Nitrogen Abundance Determination 

We derive the N/O abundance ratio from the 
[O II]A3727/[N II]A6584 ratio and assume N/O = 
N+/0+ (Pcimbert & Costero 1969). Nava et al. (2006) 
have investigated the validity of this assumption. They 
concluded that although it could be improved upon with 
modern photoionization models, it is valid to within 
about 10%. Thus, we employ this assumption, mostly 
for the purposes of direct comparison with other studies 
in the literature. 

The nitrogen to oxygen relative abundance determina- 
tions are given in Table 2. For Leo P, log(N/0) = —1.40 
± 0.09 from the KPNO 4-m spectrum and -1.33±0.05 
from the LBT/MODS spectrum. Again, there is excel- 
lent agreement between the two independent measure- 
ments. We adopt an error weighted mean of —1.36 ± 
0.04. 

The value of N/O in Leo P is high relative to most 
determinations of N/O at very low O/H, and significantly 
higher than the very narrow plateau at —1.60 ± 0.02 
which Izotov & Thuan (1999) highlighted in their study 
of very low metallicity emission line galaxies. This result 
and its implications for nitrogen nucleosynthesis will be 
discussed in Section 5.3. 

4.3. Neon, Sulfur, and Argon Abundances 

To estimate the neon abundance, we assume that Ne/O 
= Ne++/0 ++ (Peimbert & Costero 1969). The neon to 
oxygen relative abundance determinations are given in 
Table 2. We derive log(Ne/0) = -0.72 ± 0.05 from 
the KPNO 4-m spectrum and -0.78 ± 0.04 from the 
LBT/MODS spectrum. These determinations are nearly 
identical and we adopt —0.76 ± 0.03 as our final deter- 
mination. 

To determine the sulfur and argon abundances, for di- 
rect comparison, we adopt the ionization correction fac- 



tors (ICF) of Thuan et al. (1995) from H II region pho- 
toionization models to correct for the unobserved S +3 , 
Ar +2 , and Ar +4 states. In the KPNO 4-m spectrum, no 
[S III] emission lines were detected, so we do not cal- 
culate a sulfur abundance as the ICF becomes too un- 
certain (Garnett 1989), and Ar+ 3 A7136 was beyond the 
usable wavelength coverage of the spectrum. Thus, we 
only have S/O and Ar/O relative abundance measure- 
ments from the LBT/MODS spectrum. For log(S/0) we 
obtain -1.49 ± 0.07 and for log(Ar/0) we obtain -2.00 
± 0.09. 

4.4. The He Abundance of Leo P from LBT/MODS 

The LBT/MODS spectrum of Leo P shows the detec- 
tion of several He I emission lines at high significance. 
Thus, we can use the methodology developed in Olive & 
Skillman (2004) and Aver et al. (2010, 2011) to derive a 
helium abundance. Specifically, we use a Markov Chain 
Monte Carlo (MCMC) method to efficiently explore the 
parameter space in physical conditions (i.e., temperature, 
density, neutral hydrogen fraction) and other observable 
effects on the spectrum (i.e., reddening, absorption un- 
derlying the H and He emission lines, optical depth in 
the He emission lines, collisional excitation of H and He 
emission lines). We use the electron temperature de- 
rived from the [O III] emission lines as a prior, in a 
very conservative manner (see discussion in Aver et al. 
2011), producing negligible bias and effectively eliminat- 
ing non-physical false minima. In this way we determine 
the helium abundance, the physical parameters, and the 
uncertainties derived from observations of the nebula. 

The LBT/MODS produces spectra with advantages 
over the spectra typically used to determine nebular He 
abundances. The vast majority of nebular He abundance 
determinations in the literature have been derived with a 
spectral resolution of roughly 6 - 8 A. The higher resolu- 
tions of the LBT/MODS spectrographs (2.4 and 3.4 A) 
allow us to directly measure the absorption underlying 
the H emission lines and, because the resolution is a good 
match to the intrinsic width of the emission lines, provide 
optimal sensitivity for the weak emission lines. Addition- 
ally, because LBT/MODS is a double spectrograph, all 
emission lines are observed simultaneously, unlike some 
spectrographs where one needs to make multiple obser- 
vations in order to cover the wavelength range at the ap- 
propriate spectral resolution. Table 3 presents the fluxes, 
EWs, and uncertainties for the H and He emission lines 
used in this analysis. 

The methodology developed in Aver et al. (2010, 2011) 
made use of helium emissivities calculated by Porter et al. 
(2005, 2007). Thus, in all regards, we follow the method- 
ology of Aver et al. (2010, 2011). For example, H8 and 
He I A 3889 are deblended in a self-consistent way (Aver 
et al. 2010) accounting for underlying absorption and the 
equivalent widths of the stellar absorption underlying the 
He I emission lines are assumed to be equal (Aver et al. 
2010). Recently, Porter et al. (2012) have produced a 
new set of helium emissivities, which, in principle, repre- 
sent an improvement on the older emissivities. Unfortu- 
nately, an error has been discovered in the new emissiv- 
ities (Porter et al. 2013). We are currently working on 
a comprehensive paper to determine the effect of these 
corrected new emissivities on helium abundance calcu- 
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lations, and an updated analysis of Leo P will appear 
in that paper. By using the Porter et al. (2005, 2007) 
emissivities in the present analysis, we can compare the 
derived He abundance directly to those in the literature 
in Section 5.4. 

Table 4 shows the results of two different helium abun- 
dance determinations. The first solution in column 2 
is the "standard" analysis, following Aver et al. (2011, 
2012), based on 4 H I emission lines (US, H7, H/3, and 
Ha) and 6 He I emission lines (AA 3889, 4026, 4471, 5876, 
6678, 7065). (Note that in Table 4, this is referred to as 
nine emission lines because all are referenced to H/3.) In 
column 3, we give the solution after adding in He I A5015. 
The A5015 emission line has typically not been used be- 
cause of its close proximity to [O III] A5007 results in a 
blend at lower resolution. Note, however, that the res- 
olution of the SDSS spectra allowed Izotov et al. (2007) 
to measure the He I A5015 in their very large sample of 
low metallicity H II regions. 

Table 4 shows stable, well-constrained solutions for the 
physical parameters and the He abundance. As expected, 
and consistent with previous He abundance analyses, the 
derived electron temperature of 17,100 K is in excellent 
agreement with the electron temperature derived from 
the [O III] emission lines (17,350 K), and significantly 
higher than the estimated low ionization zone tempera- 
ture (14,530 K). The electron density is in the low den- 
sity limit consistent with the [S II] and [O II] emission 
lines. The reddening, C(H/3) is in excellent agreement 
with that derived from the H lines alone (for this anal- 
ysis we use the raw fluxes, not corrected for reddening, 
so the reddening determinations are independent). The 
solution for the underlying H absorption is in excellent 
agreement with the number previously derived from the 
H lines alone. Because the favored value for the optical 
depth in the He lines is low, the solution is free from as- 
sumptions about how the optical depth effects are mod- 
eled. Also, the favored solution shows little evidence for 
neutral H, so the effects of collisional excitation of the 
lower Balmer lines are also negligible. The \ 2 value of 3.3 
(confidence level of 93%) for the "standard" analysis is 
well within the cut-off of 4.0 per degree of freedom (con- 
fidence level of 94.5%) for goodness of fit recommended 
by Aver et al. (2012). 

With the addition of the A5015 He I emission line, the 
solution for the He abundance is essentially unchanged 
(as might be expected since the line has an observational 
uncertainty of 20%). However, there is a small decrease 
in the EW of the underlying He absorption. The ad- 
dition of A5015 is expected to have the greatest impact 
on the EW of the underlying He absorption because it 
is an intrinsically weak He I line. Thus, like A4026, it 
is very sensitive to underlying absorption; however, be- 
cause A5015 is a singlet line (unlike A4026, which is a 
triplet line), it is not sensitive to radiative transfer ab- 
sorption. Note that Izotov et al. (2007) have compared 
the A5015 emission line strength to A6678 in their very 
large sample of low metallicity H II regions, and find 
generally lower ratios of A5015/A6678 when compared to 
the theoretically expected values. They interpret this 
as a departure from the standard assumption of case B 
(resonance lines are optically thick). For this quality of 
spectrum, (without detection of He I A3864) it is not pos- 
sible to further test the possibility of a departure from 



the case B assumption. Our detection of A7281, though 
not strong, is in good agreement with the case B theo- 
retical prediction. Note that if A5015 were weaker due to 
the departure from case B, that we would expect an in- 
crease in the detected underlying absorption as opposed 
to the slight decrease which results. Thus, we choose to 
present an analysis including A5015. Adding A5015 re- 
sults in a slightly higher \ 2 of 3.5, and a significantly 
lower \ 2 per degree of freedom of 1.75 (corresponding to 
a 83% confidence level). 

Two possible concerns for converting a He + /H + ra- 
tio into a He/H ratio are the presence of Hc + ^ or the 
presence of neutral He. The He II A4686 emission line 
is not detected in our spectra at high confidence, indi- 
cating that any He ++ contribution is negligible. Pagel 
et al. (1992) promoted the use of the "radiation softness 
parameter", n = (0+/S+)(S ++ /0 ++ ) (Vilchez & Pagel 
1988), as an indication whether a correction for neutral 
He is warranted. Based on photoionization models with 
model stellar atmospheres, Pagel et al. (1992) concluded 
that for values of log (17) < 0.9 the correction of neutral 
He was negligible. From the LBT/MODS spectrum, r] 
= 1.84 ± 0.64, which corresponds to log(rj) — 0.26, well 
below the region where one needs to make corrections for 
neutral He. Thus, we will adopt He/H = He+/H+. The 
value of Hc+/H+ = 0.0837^;^ from the last column 

of Table 4 converts to a value of 0.2509±g;{£^ for the He 
mass fraction. 
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Fig. 4. — The newly derived oxygen abundance and estimated 
absolute magnitude for Leo P (from Rhode et al. 2013, estimating 
a distance of 1.75 ± 0.25 Mpc) compared to the B-band luminos- 
ity/oxygen abundance relationship for dwarf star forming galaxies 
from Berg et al. (2012). For the Berg et al. sample, only galaxies 
with reliable distances (e.g., TRGB distances) and direct abun- 
dances were included (i.e., the "Combined Select" sample). The 
solid line is the regression from Berg et al. (2012). 

5. THE SIGNIFICANCE OF THE CHEMICAL 
ABUNDANCES IN LEO P 

5.1. The Luminosity- Metallicity (L-Z) Relationship and 

Leo P 

Although only a single H II region is observed in Leo P, 
we will assume that the derived abundances are represen- 
tative of the ISM in the galaxy. This follows from the 
observation that metallicity gradients and variations are 
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observed to be small or non-existent in low-mass galax- 
ies (e.g., Skillman et al. 1989a; Kobulnicky & Skillman 
1996, 1997; Lee et al. 2006b; Croxall et al. 2009; Berg et 
al. 2012). 

Berg et al. (2012) have compiled a sample of nearby 
dwarf galaxies with accurate distances measured from ei- 
ther the tip of the red giant branch (TRGB) or Cepheid 
variable stars and oxygen abundances measured from 
the "direct" method (their "Combined Select" sample). 
They have shown that there is relatively low scatter in 
the relationships between the oxygen abundance with the 
B-band luminosities, the 4.5 ^.m luminosities, and the 
masses inferred from infrared luminosity and color. The 
L-Z relationship is best understood as a fundamental re- 
lationship between galaxy mass and chemical abundance 
(Tremonti et al. 2004). The low abundances in dwarf 
galaxies reflect inefficient chemical evolution due to large 
gas mass fractions and the reduced capacity to retain 
newly synthesized heavy elements. 

In Figure 4, using the estimated B-band luminosity of 
Leo P from Rhode et al. (2013), we compare the posi- 
tion of Leo P with the B-band luminosity/oxygen abun- 
dance relationship determined from the "Combined Se- 
lect" sample by Berg et al. (2012). The uncertainty in the 
B-band luminosity of Leo P is dominated by the uncer- 
tainty in its distance. Rhode et al. (2013) have estimated 
a distance of roughly between 1.5 and 2.0 Mpc from ob- 
servations of the resolved stars. This is supported by the 
distance estimate of 1.3j^'5 Mpc derived by Giovanelli 
et al. (2013) from the baryonic Tully-Fisher relationship. 
Adopting a distance estimate of 1.75 ± 0.25 Mpc and the 
photometry from Rhode et al. (2013) results in a B-band 
luminosity of -8.94 ± 0.30 for Leo P. 

Figure 4 shows the excellent agreement between the 
measurements for Leo P and the L-Z relationship of Berg 
et al. (2012). In this regard, Leo P represents an ex- 
tension of this relationship to very low luminosities for 
actively star forming dwarf galaxies. It would be very 
valuable to have a more secure distance to Leo P in or- 
der to confirm this result. Recent ground-based imaging 
observations (K. McQuinn, private communication) have 
been obtained that reach substantially deeper than the 
WIYN images from Rhode et al. (2013). The RGB is 
definitively detected in these images, and a preliminary 
analysis of these new data produce a distance determi- 
nation that is consistent with the range quoted above. A 
thorough analysis and presentation of these new obser- 
vations will be forthcoming. In fact, future observations 
of distances to other low metallicity objects would be 
equally valuable 11 . 

The Berg et al. (2012) sample, based on the nearest 
galaxies, is essentially a volume limited sample. Thus, 
the relationships from this sample should be represen- 
tative of typical galaxies (with the exception of galaxies 
experiencing severe environmental influences). The im- 
plication is that the ISM oxygen abundance of a star 

11 For example, DDO 68 (also known as UGC 5340), with an 
estimated distance of ~ 6 Mpc (from brightest stars, Makarova & 
Karachcntscv 1998) and 10 Mpc (after correcting the velocity for 
the Local Void, Pustilnik & Tcpliakova 2011), does not yet have 
a published TRGB distance. Note that Hubble Space Telescope 
observations have been acquired (program HST-GO-11578) and a 
preliminary distance of 12.1 ± 0.7 Mpc has been derived (A. Aloisi, 
private communication) . 



forming galaxy is a very strong function of the stellar 
mass, down to the least massive star forming galaxies 
known. Given that the luminosity function for galaxies 
predicts large numbers of these very low mass galaxies, 
very low metallicity galaxies like Leo P should be quite 
common. However, the requirement that they possess 
relatively high surface brightness H II regions for abun- 
dance analysis may limit their detectability significantly. 

5.2. The a Abundances and Leo P 

Figure 5 shows the Ne/O, S/O, and Ar/O relative 
abundances in Leo P as compared with the relative abun- 
dances in emission line galaxies as studied by Izotov & 
Thuan (1999). For clarity, data points with uncertainties 
of greater than 0.1 dex have been excluded. Because a 
element production is thought to be dominated by core 
collapse supernova production, under the assumption of 
a universal mass function (IMF), the ratios of the al- 
pha elements are expected to be constant as a function 
of metallicity. The Ne/O ratio in Leo P is consistent 
with the mean value and scatter of the very metal poor 
emission line galaxies shown in Figure 5. The S/O ra- 
tio in Leo P is consistent with the scatter seen in the 
metal poor emission line galaxies. The Ar/O for Leo P 
appears a bit higher than the typical emission line galax- 
ies, although consistent with the scatter. Given that the 
ionization correction for Ar/O (Izotov & Thuan 1999) 
is the most uncertain (the errorbars in the observations 
do not account for uncertainties in the ionization cor- 
rections), at this time we cannot conclude that there is 
strong evidence for an anomalous Ar/O ratio in Leo P. 
Note that Stevenson et al. (1993) have suggested that 
comparing the [Ar III] A7136 emission line to the [S III] 
A9069 emission line as a measurement of Ar/S may be 
a more reliable method for testing the stability of rela- 
tive Ar abundances. Our Ar/S measurement for Leo P 
of ^0.3 is elevated relative to the typical value of 0.2 
(Stevenson et al. 1993; Izotov & Thuan 1999), as indi- 
cated by Figure 5. 

It is well established that both emission line dwarf 
galaxies and relatively quiescent (or low star formation 
rate) dwarf galaxies show relatively constant values of 
Ne/O, S/O, and Ar/O (Thuan et al. 1995; van Zee et al. 
1997a; Izotov & Thuan 1999; van Zee & Haynes 2006). 
This has generally been regarded as expected under the 
assumption of a universal IMF; however, recent work re- 
vealing a trend of lower Ha-to-UV flux ratio with de- 
creasing galaxy luminosity has brought that assumption 
into question for dwarf galaxies (Hoversten & Glazebrook 
2008; Meurer et al. 2009; Lee et al. 2009; Boselli et al. 
2009). The constancy of the elemental abundance ra- 
tios has not played a significant role in this debate, but 
clearly the observed trends favor a universal IMF. Thus, 
alternative explanations of the Ha-to-UV flux ratio trend 
assuming a universal IMF (e.g., Fumagalli et al. 2011; 
Weisz et al. 2012) are likely to be favored, in agreement 
with the general conclusion of Bastian et al. (2010) that 
there is little secure evidence supporting a non-universal 
IMF. 

5.3. N/O Relative Abundances and Leo P 

Garnett (1990) first showed that the N/O ratio in low 
metallicity star forming galaxies is relatively constant as 
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Fig. 5. — The Ne/O, S/O, and Ar/O ratios in Leo P compared 
to the observed ratios in emission line galaxies from the samples 
of Izotov k. Thuan (1999) and Izotov et al. (2012). Data points 
with uncertainties of greater than 0.1 dex have been excluded for 
clarity. Note that all three panels have a y-axis range of 1 dex. 
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Fig. 6. — The N/O ratio in Leo P compared to the observed ratios 
in emission line galaxies from the samples of Izotov & Thuan (1999) 
and Izotov et al. (2012). Note that Leo P lies significantly above 
the very narrow plateau at log(N/0) Ri — 1.6 for the emission line 
galaxies with 12 + log(0/H) < 7.6 previously noted by Izotov & 
Thuan (1999). This indicates the likely presence of a secondary 
component of N in Leo P. 

a function of O/H (with a mean value of log(N/0) = 
-1.46t°;i3 for tnese "plateau" objects). Although the 
mean is well defined, Garnett (1990) also pointed out 
that the scatter in N/O at a given O/H is larger than 
can be accounted for in terms of observational errors. He 
suggested that the relatively large scatter in N/O at low 
O/H could be understood in terms of the different deliv- 
ery times of N and O. That is, O production is under- 



stood to be dominated by primary nucleosynthesis from 
massive stars and delivered early after a star formation 
event; N can be produced by massive stars (and deliv- 
ered early with the primary oxygen) and by intermediate 
mass stars (and delivered relatively later). N can have 
both primary and secondary origins. Thus, in a galaxy 
where the N/O ratio has been elevated from its initial 
N/O ratio from massive stars by contributions from de- 
layed nitrogen from intermediate mass stars, a burst of 
star formation will drive a galaxy to higher values of O 
and lower values of N/O, but later contributions from 
intermediate mass stars will raise the N/O at a constant 
O/H. At higher metallicities (i.e., 12 + log(0/H) > 8.0), 
the average N/O rises with increasing O/H, and this is 
thought to be indicative of the increasing influence of sec- 
ondary nitrogen production (Pagel 1985). This hypoth- 
esis has been supported quantitatively by simple models 
of galaxy chemical evolution by, e.g., Matteucci & Tosi 
(1985), Pilyugin (1993), and Fields & Olive (1998). 

In a subsequent study of blue compact dwarf galax- 
ies, Izotov & Thuan (1999) drew attention to a different 
plateau with a very small dispersion (0.02 dex) in log 
(N/O) (with a central value of —1.60) in the extremely 
metal-poor galaxies with 12 + log(0/H) < 7.6. They 
proposed that the absence of time-delayed secondary pro- 
duction of N (and C) is consistent with the scenario that 
extremely metal-poor galaxies are now undergoing their 
first burst of star formation, and that they are therefore 
young, with ages not exceeding 40 Myr. They further 
asserted that this countered the commonly held belief 
that C and N are produced by intermediate-mass stars 
at very low metallicities (as these stars would not have 
yet completed their evolution in these lowest metallicity 
galaxies) . 

Later, van Zee & Haynes (2006) derived an average of 
log(N/0) = —1.41 for a sample of isolated dwarf irreg- 
ular galaxies. Their sample included extremely metal- 
poor objects, but there was no evidence for the very nar- 
row plateau in N/O at 12 + log(0/H) < 7.6 observed 
in blue compact galaxies by Izotov & Thuan (1999). In- 
terestingly, all of the extremely metal poor H II regions 
had N/O values in agreement with the average value of 
log(N/0) = -1.41. Nava et al. (2006) revisited the ob- 
served N/O plateau with a large set of objects and deter- 
mined a mean value for the N/O plateau of log (N/O) = 
— 1.43 with a standard deviation of 1q Q84- They further 
concluded from a \ 2 analysis that only a small fraction 
of the observed scatter in N/O is intrinsic. Based on 
a Monte Carlo analysis of the scatter in the N/O ver- 
sus O/H diagram, Henry et al. (2006) concluded that 
one could not distinguish between a delayed nitrogen hy- 
pothesis or the hypothesis that nitrogen is produced by 
massive stars alone at low metallicity. They did con- 
clude that allowing galaxy ages below 250 Myr could not 
explain the plateau morphology. 

Figure 6 shows the N/O relative abundance in Leo P as 
compared with the N/O relative abundances in emission 
line galaxies as studied by Izotov & Thuan (1999) and a 
new sample available from Izotov et al. (2012). For clar- 
ity, the emission line galaxy samples have been trimmed 
of objects with uncertainties in log(N/0) of more than 
0.1 dex. Clearly, Leo P lies well above the very narrow 
plateau identified in extremely metal poor emission line 
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galaxies by Izotov & Thuan (1999). 
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Fig. 7. — The N/O ratio in Leo P compared to the observed ratios 
in emission line galaxies from the samples of Izotov & Thuan (1999) 
and Izotov et al. (2012) and the nearby dwarf galaxies from the 
sample of Berg et al. (2012). Note the larger spread in log(N/0) 
at low values of log(0/H) for the nearby galaxies sample and that 
the scatter in log(N/0) at higher values of log(0/H) for the nearby 
galaxies is comparable to that of the emission line galaxies. 

On the other hand, the N/O abundance for Leo P lies 
within the scatter observed by Garnett (1990) and van 
Zee & Haynes (2006). In Figure 7, we have added the 
galaxies from the sample of Berg et al. (2012) (again 
limited to those with uncertainties in log(N/0) less than 
0.1 dex). In this diagram, the evidence for a very nar- 
row plateau at values of 12 + log(0/H) < 7.6 is not 
seen. Thus, the very narrow plateau seen in the emis- 
sion line galaxies cannot be simply a function of metal- 
licity. van Zee et al. (2006) and Berg et al. (2012) point 
out that there are nearby galaxies with values of 12 + 
log(0/H) < 7.6 and detailed star formation histories (de- 
rived from Hubble Space Telescope observations of their 
resolved stars) that clearly show that the bulk of their 
star formation occurred well before the last 40 Myr (i.e., 
Leo A and GR 8, Tolstoy et al. 1998; Cole et al. 2007; 
Dohm-Palmer et al. 1998; Weisz et al. 2011). This implies 
that the young galaxy hypothesis is not a valid explana- 
tion for the plateau in N/O at low metallicity. 

With the addition of Leo P, it now looks as though 
even the most metal poor galaxies show a range in N/O. 
If the scatter in the log(N/0) vs 12 + log(0/H) diagram 
is due to the time delay between producing oxygen and 
secondary nitrogen as proposed by Garnett (1990), then 
it could be that the N/O in Leo P is offset to a higher 
value due to a long period of relative quiescence in which 
the nitrogen production is allowed to essentially complete 
before the next round of oxygen production. That two of 
the three very low metallicity nearby galaxies in Figure 7 
lie above the very narrow plateau is interesting. Note 
that under the interpretation that the scatter is due to 
secondary nitrogen production, then the lowest values of 
N/O (not the mean value) represent the ratio of N/O in 
primary nucleosynthesis by massive stars. 

van Zee & Haynes (2006) looked at several variables for 
their possible influence on N/O abundance. In particu- 
lar, they found a correlation between N/O and color, in 
the sense that redder galaxies have higher N/O. This sup- 



ports the delayed release hypothesis of Garnett (1990), as 
currently quiescent galaxies would be expected to have, 
on average, higher values of N/O. Berg et al. (2012) re- 
covered this trend in their sample. For the B-V color 
of 0.36 for Leo P (Rhode et al. 2013), the relationship 
between B-V color and N/O determined by Berg et al. 
(2012) predicts a value of log(N/0) = -1.50. Although, 
formally, this is not consistent with our value of —1.36 ± 
0.04, Leo P falls well within the scatter in the relation- 
ship, which has an intrinsic dispersion of 0.14 dex, and 
thus, follows the trend. 

The question raised by the N/O versus O/H diagnostic 
diagram is whether emission line galaxies and relatively 
quiescent dwarf galaxies occupy different sections of this 
diagram and therefore represent different chemical evo- 
lution paths. Within the literature there are proponents 
of the view that emission line galaxies (or blue com- 
pact galaxies) are intrinsically different from relatively 
quiescent galaxies, and there are others who allow for 
the interpretation that emission line galaxies represent a 
phase that many, or perhaps all, dwarf galaxies can pass 
through. In the former interpretation, the comparison 
of relatively quiescent dwarfs with emission line galax- 
ies could be viewed as comparing two different families 
of galaxies. The emission line galaxies could represent 
a more homogeneous sample due to the fully populated 
IMFs normally associated with massive star formation 
events. For the latter interpretation, all of the points in 
the diagnostic diagram are directly comparable, although 
the lower SFRs of relatively quiescent dwarfs are often 
associated with under-abundances of massive stars due 
to stochastic effects. Unfortunately there are few high- 
quality spectra available for extremely low-metallicity 
galaxies. Adding new observations remains a worthwhile 
enterprise. The cause of the very narrow plateau in N/O 
in very low metallicity emission line galaxies is still lack- 
ing a definitive explanation, and we will return to this in 
Section 6.2. 



5.4. Leo P and the Primordial Helium Abundance 

Next to the cosmic microwave background radiation, 
standard big bang nucleosynthesis (SBBN) is the most 
robust probe of the early universe available (e.g., Walker 
et al. 1991). Furthermore, using the precise baryon den- 
sity as determined by WMAP (Komatsu et al. 2011), 
SBBN has effectively become a parameter-free theory 
(Cyburt et al. 2002). As such, one can use SBBN to 
make relatively precise predictions of the initial light el- 
ement abundances of D, 3 He, 4 He, and 7 Li. Therefore, 
an observational determination of these abundances be- 
comes a test of the concordance between SBBN theory 
and the analyses of microwave background anisotropies. 
To test these predictions, the observed abundances must 
be determined with high precision. Unfortunately, there 
is a logarithmic relationship between the baryon to pho- 
ton ratio, rj, and the primordial helium abundance, Y p . 
Thus, any meaningful test of the theory requires a deter- 
mination of Y p to an accuracy of a few percent. The 7- 
year WMAP value for 7] is (6.16±0.15) x 10~ 10 (Komatsu 
et al. 2011). Using the procedures of Cyburt et al. (2008), 
and assuming a neutron mean life of 885.7 ±0.8 s (Naka- 
mura et al. 2010), this translates to Y p = 0.2483±0.0002, 
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Fig. 8. — The helium mass fraction (Y) and oxygen abundance for 
Leo P compared to the abundances in emission line galaxies from 
the sample of Izotov et al. (2007) as analyzed by Aver et al. (2012, 
Averl2) and the sample from Peimbert et al. (2007, PLP07). The 
single line is the regression to the data from Aver et al. (2012). Note 
that the low value of O/H and the comparable error in Y for Leo P 
make this an important contribution in the determination of the 
primordial He abundance. Two of the galaxies are common to both 
comparison samples, I Zw 18 (at low O/H) and Haro 29 (also known 
as I Zw 36 and Mrk 209, at intermediate O/H) and their points 
are connected by dotted lines. The narrow band marked WMAP7 
is the range of values (± la) estimated for the primordial helium 
abundance following the calculation by Cyburt et al. (2008) from 
the 7-year WMAP value for the baryon-to-photon ratio (Komatsu 
et al. 2011) and assuming the neutron mean life from the Particle 
Data Group collaboration (Nakamura et al. 2010). 

a relative uncertainty of only 0.08% 12 . 

To date, the favored independent method of determin- 
ing Y p uses observations of low metallicity H II regions in 
dwarf galaxies. By fitting the helium abundance versus 
metallicity, one can extrapolate back to very low metal- 
licity, corresponding to the primordial helium abundance 
(Peimbert & Torres- Peimbert 1974). The oxygen to hy- 
drogen ratio, O/H, commonly serves as a proxy for metal- 
licity. The difficulties in calculating an accurate and pre- 
cise measure of the primordial helium abundance are well 
established (Olive & Skillman 2001, 2004; Aver et al. 
2011; Izotov et al. 2007; Peimbert et al. 2007). 

Figure 8 shows the helium mass fraction (Y) and oxy- 
gen abundance for Leo P compared to the abundances 
in emission line galaxies from the sample of Izotov et al. 
(2007) as analyzed by Aver et al. (2012) and the sample 
of very high quality spectra analyzed by Peimbert et al. 
(2007). The two comparison data sets should be directly 
comparable in that they use identical atomic data for 
converting relative emission line fluxes into abundances 
(e.g., Porter et al. 2005, 2007). However, the methodolo- 
gies for determining physical conditions vary somewhat. 
For example, two of the galaxies are common to both 
comparison samples, I Zw 18 and Haro 29 (also known 
as I Zw 36 or Mrk 209). These two galaxies are connected 
by solid lines in Figure 8, and one can see that the val- 
ues for O/H are significantly different (due primarily to 
assumptions about the effects of temperature variations 
within the nebulae). We are not advocating that the dif- 

12 Note that the Planck Collaboration has just published a 
new prediction of the primordial Helium abundance of 0.24771 
± 0.00014, corresponding to an uncertainty of 0.06% (Ade et al. 
2013). 



ferent samples be combined, but are only conducting the 
comparison to show the potential impact of Leo P on 
determinations of Y p . 

Note that the low value of O/H and uncertainties in 
He/H which are comparable to the other best studied 
nebulae in the literature make Leo P an important con- 
tribution in the determination of Y p . To demonstrate 
this, we calculated a regression of Y on O/H by adding 
our Leo P He abundance determination to the sample 
from Aver et al. (2012). This reduced the intercept (Y p ) 
by roughly 0.3% from 0.2534 to 0.2527 and the error on 
the intercept by roughly 8% from 0.0083 to 0.0076. The 
resulting regression of 

Y = 0.2527 ±0.0076 + 61 ±92(0/ H) (5) 

has a x 2 of 2.9. 

The value of Y p implied by the 7-year WMAP obser- 
vations is indicated in Figure 8. Note, especially, that 
our determination of Y for Leo P is in excellent agree- 
ment with the two determinations for Y in I Zw 18 and 
that all three are in excellent agreement with the 7-year 
WMAP prediction for Y p . Recently, Izotov & Thuan 

(2010) reported a significant difference between the value 
of Y p derived from observations of metal poor H II re- 
gions and the value calculated from the cosmic microwave 
background observations and suggested this was evidence 
for non-standard Big Bang Nucleosynthesis. From the 
present analysis of the most metal poor H II regions, we 
see no motivation for non-standard Big Bang Nucleosyn- 
thesis. 

6. LEO P AND XMD GALAXIES 
6.1. The Search for Extremely Low Metallicity Galaxies 

In an overview of the most metal-poor galaxies, Kunth 
& Ostlin (2000) defined "very metal deficient" galaxies 
as those with metallicities of ten percent of the solar 
value or less. At the time, the solar value for the oxygen 
abundance was generally accepted to be 12 + log(0/H) 
= 8.91. In the last decade, there has been a consensus 
that the solar oxygen abundance is lower than previously 
thought (e.g., 12 + log(0/H) = 8.69 ± 0.05, Asplund et 
al. 2009) and recent papers have adopted the definition of 
an "extremely metal deficient" (XMD) galaxy as having 
12 + log(0/H) < 7.65 (e.g., Kniazev et al. 2003; Pustilnik 
& Martin 2007; Kakazu et al. 2007; Ekta et al. 2008; 
Brown et al. 2008; Ekta & Chengalur 2010b). 

Also, in the last decade, there have been several pro- 
grams with the aim of discovering more XMD galaxies. 
The majority of these programs have concentrated on 
surveying emission line galaxies. Examples include Kni- 
azev et al. (2000); Ugryumov et al. (2003); Kniazev et al. 
(2003); Melbourne et al. (2004); Izotov et al. (2006); Pa- 
paderos et al. (2008); Brown et al. (2008); Guseva et al. 

(2011) . I Zw 18 (Searle & Sargent 1972) and SBS 0335- 
052 (Izotov et al. 1990) represent prototypes for XMD 
galaxies that can be found in emission line surveys. Re- 
cently, direct oxygen abundance measurements have been 
extended to higher redshifts (Hu et al. 2009; Xia et al. 
2012). It is interesting that while recent surveys have 
significantly enlarged the number of metal-poor galax- 
ies, only a few XMD galaxies with oxygen abundances 
below ~5% of the solar value are known, and no galax- 
ies with abundance below 12 + log(0/H) ss 7 have been 
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discovered in the Local Volume. 

Despite the many surveys of emission line galaxies, the 
number of XMD galaxies remains small. Searching via 
emission line surveys provides an extremely low yield. As 
reported in Izotov et al. (2012), of one million SDSS spec- 
tra, 13,000 emission line objects have detectable [O III] 
A4363, and, of these, there are only 15 candidates with 
12 + log(0/H) < 7.35. 

An alternate approach to finding low metallicity galax- 
ies relies on the fundamental L-Z relationship between 
galaxian stellar luminosity (or mass) and abundance for 
low redshift dwarf galaxies as discussed in Section 5.1 
(e.g., Peimbert & Spinrad 1970; Lequeux et al. 1979; 
Skillman et al. 1989a; Lee et al. 2006a; Berg et al. 2012). 
Because of the galaxy luminosity function, XMD galax- 
ies are very numerous. However, in part because of their 
small sizes, XMD galaxies with luminous, high surface 
brightness star forming regions are rare. The result is 
an extreme paucity of XMD galaxies in emission line 
surveys. Many XMD galaxies have been found in sur- 
veys of nearby, low luminosity galaxies. XMD galax- 
ies such as Leo A (Skillman et al. 1989a; van Zee et 
al. 2006), SagDIG (Skillman et al. 1989b; Saviane et al. 
2002), UGCA 292 (van Zee 2000), DDO 68 (Pustilnik 
et al. 2005), and two XMD galaxies in the Lynx-Cancer 
void (Pustilnik et al. 2011) have been found in this way. 
Spectroscopic follow-up of Ha surveys of low luminosity 
galaxies and galaxies discovered in blind H I surveys con- 
tinue to be very promising in this regard (e.g., Cannon et 
al. 2011). Surveys like ALFALFA have discovered hun- 
dreds of objects whose properties may be confirmed by 
subsequent follow-up observations to be consistent with 
those of nearby, low-luminosity dwarf galaxies. A ded- 
icated program to measure the nebular abundances of 
these galaxies may be very fruitful in discovering new 
XMD galaxies. 
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M(B) 

Fig. 9. — Figure 4 is reproduced showing the positions of the 
famous XMD blue compact galaxies I Zw 18 and SBS 0335-052 
relative to the B-band luminosity /oxygen abundance relationship 
for dwarf star forming galaxies from Berg et al. (2012). The 
XMD galaxies with disturbed kinematics as observed by Ekta et al. 
(2006, 2008); Ekta & Chengalur (2010a) are also shown. Note the 
large offsets of the XMD emission line galaxies from the luminos- 
ity/oxygen abundance relationship defined by dwarf star forming 
galaxies. 



6.2. XMD Galaxies and the Infall of Metal-Poor Gas 

It is well-established that many of the XMD emission 
line galaxies lie well off the L-Z relationship for dwarf 
galaxies. This includes the two best-known XMD galax- 
ies, I Zw 18 and SBS 0335-052. In Figure 9 we have re- 
produced the L-Z diagram shown in Figure 4, and added 
I Zw 18 and SBS 0335-052. For I Zw 18, we have used 
the new distance of Aloisi et al. (2007) and the abun- 
dances for the NW component from Skillman & Kenni- 
cutt (1993) which is identical to that reported by Izotov 
et al. (1999). For SBS 0335-052, we have plotted the 
eastern (luminous) component assuming a distance of 54 
Mpc and used the oxygen abundance from Izotov et al. 
(1999). In Figure 9 it is clear that these two galaxies lie 
off the L-Z relationship by about 7 magnitudes (or 0.7 
dex in log(0/H)). These offsets are much too large to 
be attributable to increases in the luminosity due to the 
present starburst. 

Based on new H I observations of XMD galaxies, Ekta 
& Chengalur (2010b) have proposed that inflow of metal- 
poor gas from the outskirts of a galaxy to the central 
star- forming regions (in the case of interaction), or cold 
gas accretion is the probable cause for the observed low 
emission-line metallicities of XMD galaxies. This hy- 
pothesis is very appealing because it solves two prob- 
lems simultaneously. It explains how small systems can 
achieve such high star formation rates that they become 
visible as emission line galaxies. It also explains why 
so many of these galaxies do not follow the metallic- 
ity/luminosity relationship observed for normal dwarf 
galaxies in volume limited surveys. 

Note that it is also possible to have very low abun- 
dances in galaxies through metal-enriched outflows. This 
is not as appealing as the Ekta & Chengalur (2010b) hy- 
pothesis because this would require that the majority of 
nucleosynthesis in the XMD emission line galaxies is as- 
sociated with strong bursts of star formation resulting in 
galactic winds. Furthermore, this does not explain why 
other galaxies at similar masses have been able to enrich 
to higher metallicities. 

In Figure 9, we have also added in the XMD galax- 
ies with disturbed kinematics as observed by Ekta et 
al. (2006, 2008) and Ekta & Chengalur (2010a). Here 
one can see similarly large offsets of the XMD galax- 
ies from the luminosity/oxygen abundance relationship 
defined by dwarf star forming galaxies. Although off- 
sets from the L-Z relationship for XMD emission line 
galaxies are well known, what Ekta and collaborators 
have added is the observation that all of these galaxies 
have disturbed kinematics indicative of either interac- 
tion or infall. In support of this hypothesis, I Zw 18 and 
SBS 0335-052 are both known to be associated with dis- 
turbed H I kinematics (van Zee et al. 1998; Lelli et al. 
2012; Pustil'nik et al. 1997; Ekta et al. 2009). 

It would be very satisfying if this hypothesis could also 
explain the very narrow plateau seen in the N/O abun- 
dances in XMD emission line galaxies. One could hy- 
pothesize that if the XMD galaxies are being polluted 
by relatively pristine gas, and if that pristine gas had 
a single value of N/O, then XMD emission line galaxies 
would tend to have similar values of N/O. Unfortunately, 
absorption line studies of damped Lyman-a systems do 
not show a universal value of N/O; in fact they show a 
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very large range in N/0 at low values of O/H (Pcttini et 
al. 2008) (which, interestingly, has been used as evidence 
for a non-universal IMF, Tsujimoto & Bekki 2011). If all 
emission line XMD galaxies are the products of strong 
interactions and infall, then it makes the very narrow 
plateau seen in the N/O abundances in XMD emission 
line galaxies even more difficult to understand. 

The discovery of Leo P supports the hypothesis that 
XMD galaxies are of two different types. XMD emis- 
sion line galaxies represent starbursts which are likely 
triggered by interaction with a companion or infall from 
an H I cloud, while relatively "quiescent" (i.e., galaxies 
with low, but non-zero star formation rates) dwarf XMD 
galaxies are formed through normal evolution of very low 
mass galaxies. 

7. CONCLUSIONS 

We have presented KPNO 4-m and LBT/MODS spec- 
troscopic observations of an H II region in the nearby 
dwarf irregular galaxy Leo P which was discovered re- 
cently in the Arecibo ALFALFA survey. The results 
from both observations are in excellent agreement. The 
KPNO and LBT observations yield accurate measure- 
ments of the temperature sensitive [O III] A4363 line and 
a "direct" oxygen abundance of 12 + log(0/H) = 7.17 ± 
0.04. This oxygen abundance is among the lowest ever 
measured in an H II region. Thus, Leo P is an extremely 
metal deficient (XMD) galaxy. 

For its estimated luminosity, Leo P is consistent with 
the relationship between luminosity and oxygen abun- 
dance seen in nearby dwarf galaxies. Future observations 
of distances of nearby XMD galaxies will help to con- 
firm that the relationship between luminosity and oxy- 
gen abundance extends down to the least massive star 
forming galaxies known. 

Leo P shows normal a element abundances (Ne/O, 
S/O, and Ar/O) when compared to other XMD galax- 
ies. These well-defined trends in the lowest metallicity 
galaxies are supportive of the hypothesis of a universal 
initial mass function. 

The N/O ratio in Leo P is elevated relative to sim- 
ilarly low metallicity emission line galaxies, but within 
the scatter seen in observations of large samples of dwarf 
galaxies. The elevated N/O and optical color of Leo P 
are consistent with the "delayed release" hypothesis for 
N/O abundances. 

The high signal/noise ratio and spectral resolution of 
the LBT/MODS spectrum allow us to derive a helium 
abundance with a precision of ~ 5% . We derive a helium 
mass fraction of 0.2509t£; 

0123 which compares well with 
the WMAP + BBN prediction of 0.2483 ± 0.0002 for 
the primordial helium abundance. Supplementing the 
Leo P observations with similar quality spectra of other 
XMD galaxies will allow a future determination of the 
primordial helium abundance with an accuracy of less 
than two percent. 

The discovery of Leo P suggests that surveys of very 
low mass galaxies compete well with emission line galaxy 
surveys for finding XMD galaxies. With present and fu- 
ture blind H I surveys producing large numbers of low 
mass galaxies, it should be possible to greatly increase 
the number of studies of abundances in XMD galaxies. 

It is possible that XMD galaxies may be divided into 



two classes. The XMD emission line galaxies are rare. 
This may be attributable to the uncommon occurrence 
of starbursts triggered by infall or interaction. Because 
of the galaxy luminosity function (which indicates large 
numbers of very low mass galaxies) and the relation- 
ship between galaxy stellar mass and oxygen abundance, 
XMD galaxies with very low chemical abundances due to 
their intrinsically small masses should be common. 

Our ability to observe the XMD galaxies with intrinsi- 
cally small masses may be limited by their low star for- 
mation rates and the relatively short lives of the massive 
stars that produce the observable H II regions. Alter- 
nately, their intrinsically low luminosities may prevent 
them from inclusion in most galaxy catalogs. For exam- 
ple, the Ha imaging survey of galaxies in the local 11 
Mpc volume by Kennicutt ct al. (2008) had an appar- 
ent magnitude cut-off of tub < 15. With an mg of 17.3 
(Rhode et al. 2013), Leo P is more than 2 magnitudes 
fainter than this cut-off. Blind H I studies may be the 
only efficient way to detect these low luminosity XMD 
galaxies. 
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TABLE 1 

Emission-Line Intensities and Equivalent 
Widths for Leo P 



J(A)//(H/3) 



Ion 

[O II] A3727 

H12 A3750 

Hll A3771 

H10 A3798 

He I A3820 

H9 A3835 

[Ne III] A3868 

He I+H8 A3889 

[Ne III]+H7 A3968 

He I A4026 

H<5 A4101 

H7 A4340 

[O III] A4363 

He I A4471 

H/3 A4861 

He I A4921 

[O III] A4959 

[O III] A5007 

He I A5015 

He I A5876 

[S III] A6312 

[N II] A6548 

Ho A6563 

[N II] A6584 

He I A6678 

[S II] A6717 

[S II] A6731 

He I A7065 

[Ar III] A7136 

He I A7281 

[O II] A7320 

[O II] A7330 

P13 A8665 

P12 A8750 

Pll A8863 

P10 A9015 

[S III] A9069 

P9 A9229 

[S III] A9532 

P8 A9546 



KPNO 4-m LBT/MODS 



0.467 ± 0.038 

0.040 ± 0.010 
0.045 ± 0.010 

0.062 ± 0.009 
0.118 ± 0.011 
0.195 ± 0.015 
0.181 ± 0.014 

0.253 ± 0.015 
0.473 ± 0.020 
0.035 ± 0.007 
0.034 ± 0.005 
1.000 ± 0.022 

0.468 ± 0.010 
1.399 ± 0.031 

0.094 ± 0.006 

0.008 ± 0.004 
2.699 ± 0.197 
0.021 ± 0.005 
0.026 ± 0.006 
0.047 ± 0.007 
0.027 ± 0.006 



0.465 
0.038 
0.043 
0.055 
0.010 
0.072 



0.017 
0.007 
0.007 
0.008 
0.007 
0.007 



0.038 
0.034 
1.000 
0.009 
0.482 
1.453 
0.025 
0.098 



0.005 
0.006 
0.020 
0.006 
0.010 
0.030 
0.005 
0.004 



0.106 ± 0.008 
0.197 ± 0.009 
0.186 ± 0.009 
0.011 ± 0.007 
0.268 ± 0.009 
0.452 ± 0.011 
± 
± 
± 
± 
± 
± 
± 
± 

0.007 ± 0.003 
0.009 ± 0.003 
2.750 ± 0.095 
0.025 ± 0.003 
0.028 ± 0.003 
0.036 ± 0.003 
0.027 ± 0.003 
0.023 ± 0.002 
0.026 ± 0.002 
0.007 ± 0.002 
0.006 ± 0.002 
0.007 ± 0.002 
0.007 ± 0.002 
± 
± 

0.015 ± 0.003 
0.043 ± 0.004 
0.020 ± 0.003 
0.102 ± 0.007 
0.023 ± 0.003 



0.013 
0.011 



0.002 
0.002 



-cm) 

EW(H/3(ABS)) (A) 
F(H/3) 
EW(H/3) (A) 
EW(Ha) (A) 



0.00±0.10 
2.0±2.0 
46±0.7 
143 
927 



0.09±0.04 
1.0±2.0 
34±0.7 
183 
1222 



Note. — Emission line fluxes (measured by directly 
integrating under the line profile and only using de- 
blended Gaussian profile fits and multiple component 
fits when necessary) are relative to H/3 = 1.00 and are 
corrected for reddening. The H/3 flux is given for ref- 
erence, with units of 10 -16 erg s _1 cm . Note that 
uncertainties listed in this table reflect the statistical 
uncertainties in the flux through the slit only, and do 
not account for slit losses. 
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TABLE 2 

Ionic and Total Abundances for Leo P 







KPNO 4-m 

1V1 IN vy rt 111 


LRT/MODS 








T(0++) 

TYO+I ("inferred"! 
T(S++) (inferred) 
n e (b T j 
n e (0+) 




17150 + 20 i° 

— 1390 

i ^q^n + 1 1 on 

n +235 ^ TV ,-3 
U _q cm 


17350 +"90 

— lUoU 
1 4^0 4- 1 0^0 

□! LUOU 

1UU1U _1_ -L ± '■±1/ 

fin +200 „ m -3 
bu __qq cm 

45 +45 cm -3 








(V /ft) 




n A9 -i- n i ^ vin~ 


n 1 1 J_ n no v 1 O" 
u.^ti nz u.uy xiu 


-5 






(O++/H) 




1.08 ± 0.22 xl0~ 5 


1.09 ± 0.16 xlO- 


-5 






(O/H) 




1.49 ± 0.26 xl0~ 5 


1.50 ± 0.18 xlO" 


-5 


1.46 ± 0.14 Xl0" 


-5 


log(0/H) + 12 




7.17 ± 0.07 


7.17 ± 0.05 




7.17 ± 0.04 




(N+/O+) 




4.0 ± 0.9 xl0~ 2 


4.7 ± 0.6 xl0~ : 


! 


4.4 ± 0.4 xlO" 


2 


log(N/0) 




-1.40 ± 0.09 


-1.33 ± 0.05 




-1.36 ± 0.04 




(N+/H) 




1.7 ± 1.1 X 10 


2.0 ± 0.7 X 10 








ICF 


3.58 ± 0.19 


3.67 ± 0.13 








(N/H) 




5.9 ± 1.7 xl0~ 7 


7.1 ± 1.2 xlO" 




6.7 ± 1.0 xlO" 


T 


log(N/H) + 12 




5.77 ± 0.11 


5.85 ± 0.07 




5.82 ± 0.06 




(Ne++/0++) 




1.92 ± 0.25 xlO" 1 


1.66 ± 0.14 xlO" 


-1 


1.72 ± 0.12 xlO" 


-1 


log(Ne/0) 




-0.72 ± 0.05 


-0.78 ± 0.04 




-0.76 ± 0.03 




/■\T„ + + /TT\ 




2.0/ ± U.56 X 10 


1 on _l_ n q £ win- 
1.8U ± O.OO X 10 


-6 






ICF 


1.39 ± 0.19 


1.37 ± 0.13 








(Nc/H) 




2.87 ± 0.62 xlO- 6 


2.48 ± 0.37 xlO" 


-6 


2.58 ± 0.32 xlO" 


-6 


log(Ne/H) + 12 




6.46 ± 0.08 


6.39 ± 0.06 




6.42 ± 0.05 




(S+/H) 




0.75 ± 0.25 xl0~ 7 


0.64 ± 0.12 xlO" 


-7 






(S++/H) 






2.84 ± 0.33 xlO" 


-7 








ICF 




1.39 ± 0.12 








(S/H) 






4.84 ± 0.55 xlO" 


-7 


4.84 ± 0.55 xlO" 


-7 


log(S/H) + 12 






5.68 ± 0.47 




5.68 ± 0.47 




log(S/0) 






-1.49 ± 0.07 




-1.49 ± 0.07 




(Ar++/H) 






9.1 ± 1.3 xl0- f 


i 








ICF 




1.65 ± 0.18 








(Ar/H) 






1.50 ± 0.27 xlO" 


-7 


1.50 ± 0.27 xlO" 


-7 


log(Ar/H) + 12 






5.18 ± 0.07 




5.18 ± 0.07 




log(Ar/0) 






-2.00 ± 0.09 




-2.00 ± 0.09 





TABLE 3 

Inputs and x 2 s for. MCMC Analysis of Helium Abundance 



Emission Line 


Flux 


EW 


x 2 


+ A5015 x 2 


H<5 A4101 


0.246 ± 0.006 


26.5 ± 2.7 


0.252 


0.310 


H7 A4340 


0.431 ± 0.009 


54.6 ± 5.5 


2.751 


2.686 


H/9 A4861 


1.000 ± 0.020 


183.2 ± 18.3 






Ho A6563 


2.949 ± 0.059 


1222 ± 122 


0.191 


0.177 


He I+H8 A3889 


0.175 ± 0.007 


16.3 ± 1.6 


0.019 


0.009 


He I A4026 


0.011 ± 0.007 


1.0 ± 0.1 


0.004 


0.029 


He I A4471 


0.033 ± 0.006 


4.5 ± 0.5 


0.015 


0.034 


He I A5015 


0.026 ± 0.006 


5.0 ± 0.5 




0.211 


He I A5876 


0.103 ± 0.004 


30.2 ± 3.0 


0.001 


0.002 


He I A6678 


0.030 ± 0.003 


14.3 ± 1.4 


0.024 


0.022 


He I A7065 


0.025 ± 0.003 


12.8 ± 1.3 


0.028 


0.044 


Total x A 






3.29 


3.50 



Note. — Emission line fluxes are relative to H/3 = 1.00 and are not 
corrected for reddening (as reddening is derived in the MCMC analysis). 
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TABLE 4 

MCMC Analysis of Helium Abundances and 
Physical Conditions for Leo P 



Parameter 



Standard 



A5015 



Emission lines 
Free Parameters 
(He+/H+) 
T e (10 4 K) 
n e (cm -3 ) 
C(H/3) 
a H (EW A) 
a H e (EW A) 

T He 

n HI (1(T 4 cm" 3 ) 
X 2 

X 2 /d.o.f. 
Y 



0.0837 
1.71 



9 
8 

+0.0084 
0.0062 
+0.19 



-0.28 
1 + 20fi 

o m+ 03 

.+1.44 
*-0.94 
n +0.42 
-0.42 
1+0.66 
-0.00 
0+156 



0.94! 
0.50j 

0.00 4 



-0 

3.3 
3.3 

2509+ ' 0184 
u.zouy_ 0142 



0.0837 



10 

8 

+0.0084 
0.0054 



1.72 



+0.18 



0.27 
+201 



-1 



1-0.03 



0- IOIo.ot 

1 02+ 1 ' 43 

1 - uz -1.02 

u - 40 -0.45 
-,+0.65 
-0.00 
0+152 



o.oo: 



-0 

3.5 

1.75 

2509+ ' 0184 
u.zouy_ 0123 



